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ABSTRACT 

We present the results of the pulse phase- and luminosity-resolved spectroscopy of the 
transient X-ray pulsar V 0332-1-53, performed for the first time in a wide luminosity range 
(1 - 40) X 10^^ erg s“' during a giant outburst observed by the RXTE observatory in Dec 
2004 - Feb 2005. We characterize the spectra quantitatively and built the detailed “three- 
dimensional” picture of spectral variations with pulse phase and throughout the outburst. We 
show that all spectral parameters are strongly variable with the pulse phase, and the pattern of 
this variability significantly changes with luminosity directly reflecting the associated changes 
in the structure of emission regions and their beam patterns. Obtained results are qualitatively 
discussed in terms of the recently developed reflection model for the formation of cyclotron 
lines in the spectra of X-ray pulsars. 
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1 INTRODUCTION 


The observed emission of X-ray pulsars originates from one or two 
hot regions, located near the neutron star magnetic poles where 
the accreting matter is funneled to by the strong magnetic field 
of the neutron star. The geometry of the emission regions de¬ 
pends mostly on the accretion rate, magnetic field structure and 
stre ngth, and can be described as spots, mounds or columns (see 
e.g..lBas ko & SunvaevI 197^:lBumMd et al.|[l9^ : lBecker & Wolfl] 
Q7|; ■ 


120071: iMushtukov et alj l2014h . The geometry and the beaming 


properties of these hot regions shape the observed X-ray spectrum, 
its pulse-phase and accretion-rate dependences which can be stud¬ 
ied observationally and used to uncover intrinsic properties of the 
emission regions. 

Significant progress in X-ray astronomy has been achieved 
in the last decades with the Compton-GRO, RXTE, Bep- 
poSAX, INTEGRAL and Swift observatories, which allowed 

(see. 


e.g., reviews of 

Bildsten et al. 

1997 

; dal Fiume et al. 20001: 

Cobum et al.ll2002l; 

Filippova et al.l 

2005 

:lLutovinov & TsvgankovI 

2009al; Caballero & Wilms||2012|). A combination of excellent tim- 
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ing capabilities of the RXTE observatory with its good spectral 
characteristics allowed also to perform a pulse phase-resolved anal¬ 
ysis for a number of X-ray pulsars (see, e.g ., iKreykenbohm et ^ 
|2004 iKlochkov et al.l20o4 iFiirst et al.ll201 ll) . 

In this work we focus on the detailed pulse phase-resolved 
spectroscopy for the X-ray pulsar V 0332-1-53 in a wide dynamical 
luminosity range L 37 ^1-40 (here L 37 = Ty/lO^^ erg s~‘, as sum- 
ing the distance to the source of 7 kpc. lNegueruela et al.ll999h . It is 
important to emphasize, that such systematic analysis has not been 
carried out before. 


The X-ray pulsar V 0332-t-53 was observed by the RXTE 
observatory during a powerful outburst in 2004-2005, that pro¬ 
vided a very good coverage of different outburst stages. An ad¬ 
ditional reason to choose V 0332-t53 for the study of physical 
conditions and geometry of the emission regions in X-ray pul¬ 
sars a nd for an examinatio n of a recently developed reflection 
model dPoutanen et alj2o'l^ is a presence of the cyclotron absorp¬ 
tion line (so-called cyclotron resonant scatterin g feature, CRSF) 
and its higher harmonics in the source spectrum 1 Makishima et al.l 
Il99d : l^ttschmidt et al.ll 200 i : IKreykenbohm et"^ 2005 ). A study 
of the source pulse-averaged spectra showed that the cyclotron line 
energy, as well as parameters of the spec tral continuum, signifi- 
cantly change with the source luminosity dTsygankov et al.|[2006l : 
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Table 1. Log of RXTE observations (ordered by luminosity) used for the pulse phase-resolved analysis of the X-ray pulsar V 0332+53. 


# 

Observation 

Date 

Luminosity ^ 

Phase ^ 


ID 

MJD 

(10^^ erg s“*) 


1 

90089-11-04-04 

53360.0 

36.9 ±0.2 

m 

2 

90089-11-04-02G 

53358.6 

35.4 ±0.2 

b 

3 

90089-11-04-00G 

53356.5 

34.3 ±0.3 

b 

4 

90089-11-03-0IG 

53354.7 

32.8 ±0.3 

b 

5 

90089-11-03-03 

53352.8 

30.4 ±0.2 

b 

6 

90089-11-03-04 

53353.7 

28.5 ± 1.1 

b 

7 

90427-01-03-01 

53385.1 

22.5 ±0.8 

d 

8 

90427-01-03-05 

53387.0 

20.2 ±0.1 

d 

9 

90427-01-03-06 

53387.4 

19.6 ±0.1 

d 

10 

90427-01-03-09 

53388.4 

18.6 ±0.4 

d 

11 

90427-01-03-12 

53389.3 

18.0 ±0.5 

d 

12 

90014-01-02-13 

53390.4 

17.3 ±0.1 

d 

13 

90089-11-02-03 

53343.8 

16.3 ±0.1 

b 

14 

90089-11-02-00 

53342.8 

15.4 ±0.1 

b 

15 

90014-01-03-020 

53394.4 

15.3 ±0.1 

d 

16 

90014-01-03-03 

53395.3 

14.3 ±0.1 

d 

17 

Group 1: 
90014-01-04-00 

53398.5 

12.7 ±0.2 

d 


90014-01-04-01 

53399.6 




90014-01-04-02 

53401.4 



18 

90014-01-05-01 

53407.6 

9.7 ±0.1 

d 

19 

90427-01-04-04 

53413.8 

7.0 ± 0.2 

d 

20 

Group 2: 
90427-01-04-02 

53414.3 

6.2 ±0.1 

d 


90427-01-04-03 

53413.8 




90427-01-04-05 

53414.8 



21 

90427-01-04-01 

53416.6 

5.0 ±0.1 

d 

22 

Group 3: 
90014-01-07-01 

53419.4 

3.2 ±0.1 

d 


90014-01-07-02 

53420.7 




90014-01-07-03 

53420.7 



23 

90014-01-07-00 

53424.4 

1.7 ±0.1 

d 


" Luminosity in the 3-100 keV energy band assuming the source distance of 7 kpc 
^ Outburst phase: b - brightening, d - declining, m - near the maximum 


iMowlavi et Ii]|2006l:lTi7gankov et alJl2010l:lNakaiima et alfeOlOh 

along the outburst. This behaviour was interpreted as a result of 
changes in the height of the accretion column wher e the cyclotron 


absorption line was a s sumed to origina t e (see . e.g..lMihara 


ISchbnherr et aD l200l iNishimiiralEool I20l4 ISecker et al 


1995 


2012 , 


and references the rein). However, alterna tive scenario has recently 
been proposed bv IPoutanen et al.l l l2013l) . where the line forms in 
the illuminated atmosphere of the neutron star. In our paper, we 
focus on the latter model when discussing observed spectral varia¬ 
tions. 

According to the model calculations the cyclotron line param¬ 
eters should strongly depend on the viewin g angle, pulse phase , 
bulk velocity of the falling plasma, etc. 


Arava^IIardin^ 1994 Arava-Gochez & 


m g angle, pulse pnase , 

. fisenberg et alj|l998l: 

Harding 200d : ISerbej 


20001 : ISchbnherr et alj2007h . Parameters of the continuum should 
also depend on the viewing angle due to the significant angular 
dependence of the continuum opacit y in the strong magnetic fiel d 
near the neutron star surface (see, e.g.. lKaminker et aljl982l . ri983h . 
Thus, observations of V 0332-1-53 give us an unique possibility to 
build for the first time a detailed “three-dimensional” picture (pa¬ 
rameter vs pulse phase vs luminosity) of the source spectrum be¬ 
haviour, and ultimately uncover the physical and geometrical pic¬ 
ture of the accretion onto a magnetized neutron star. 


Note, that first results on the pulse pha se-resolved spec¬ 
trosco py for V 03324-53 were obtained by IPottschmidt et all 
( l2005h . however, available observations did not contain the PCA 
data with the required energy and time resolution, so the analysis 
was limited to the HEXTE data and, therefore, it was not trivial 
to constrain the parameters of the broadband spectrum. As a re¬ 
sult, no significant variations of spectral parameters with the pulse 
phase had been reported. Later, some preliminary results (using 
only two observations of V 03324-53) on the phase and luminos¬ 
ity depen dence of the continuum and CRS F parameters were pub¬ 
lished bv iLutovinov & Tsygankovl l l2009bh and the results on vari- 
abi lity of the iron emission line with the pulse phase and luminosity 
by iTsygankov & Lutovinotl (l2009l) . Also no attempts to interpret 
the obtained results have been made. 


2 OBSERVATIONS, SPECTRAL MODEL AND DATA 
ANALYSIS 

As already mentioned, V 03324-53 has been observed with the 
RXTE observatory durin g the so-called gi ant outburst, which 
started in the end of 2004 dSwank et al.ll2004) and lasted for about 
two months. In the current work we use the same data set as 
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Figure 1. Representative examples of the pulse phase-resolved spectra of V 0332-1-53 for two different luminosities L 37 = 36.9 (obs. 1, left) and 20.2 (obs. 
8 , right). Black symbols correspond to the spectra at phases 0.5625-0.625 (left) and 0.75-0.8125 (right), red symbols correspond to the spectra at phases 
0.9375-1.0 (left) and 0.5-0.5625 (right), the blue symbols are for the phases 0.5-0.5625 (left) and 0.375-0.4325 (right). Values of^^ and the cyclotron line 
energy are also shown for reference. Residuals from the best-fitting models (summarised in Table|3 are shown in lower panels. 


Table 2. Best-fitting parameters for three pulse phase-resolved spectra for two different luminosities of V 0332-t-53. 


Phase 

r 

•^cut 

Eo 

o-Q ro 

El 

0-1 

ri 

(d.o.f.) 



keV 

keV 

keV 

keV 

keV 







Obs.l,L37 = 36.9 





0.5625-0.625 

0.9375-1.0 

0.5-0.5625 

0.35 ± 0.07 
0.42 ± 0.05 

r» 1 0+O.O8 
^■^°-0.04 

n t+0.7 
'•■^-0.5 

7 6^^-^ 
'•^-0.4 

^•^“^-0.44 

23.72 ±0.15 
23.93 ±0.11 
23 lO^**-** 

6.57+“ ’2 1.47 ±0.08 

6 . 68 +°'« 1.93 ±0.08 

7.15+°'“ 1.38 ±0.08 

Obs. 8 ,'L 37 = 20.2 

51.6!“ 

51.6!+^ 

50-2!u 

8.2 ±3.8 

10.5 ±2.8 

1 Q "2+0.64 
^•°^-0.34 

2 56+° °^ 

262+8:^^ 

-0.30 

91.9 (90) 
97.5 (90) 
94.1 (90) 

0.75-0.8125 

0.5-0.5625 

0.375-0.4325 

0.34 ± 0.07 
0.28 ± 0.08 
0.38 ± 0.08 

0 c+0.5 

8 .6_og 

7.7 ±0.7 

25 54+*’-*^ 
26.17 ±0.27 

26 60+**'^^ 
ZO.OU Q 

6.00+“^ 1.81 ±0.07 

7.15+ll'g 1.88 ±0.13 

5.41+°'” 1.76 ±0.09 

53.4!{;| 

-0.9 

9 0+^* 

‘^•‘^-2.6 

6.4 ± 2.5 

2.37 ±0.32 

4 0+^ '* 

‘t.U-i 2 

0 '20+0.83 

111.9(107) 
90.1 (107) 
99.5 (107) 


in iTsygankov et al.l ( 1201 Ol) apart from a number of observations 
(ObsIDs. 90014-01-01-XX, 90089-11-05-XX, 90089-22-01-XX, 
90427-01-01-XX, 90427-01-02-XX), performed with time resolu¬ 
tion insufficient for the pulse-phase resolved analysis (typically 16 
s, while the source pulse period is about of 4.35 s). The used data set 
covers both brightening and declining phases of the outburst more 
or less uniformly. There is only a small gap between luminosities 
of L 37 = 22.5 and L 37 = 28.5 due to the absence of observations 
of the required timing quality. To improve the statistic for obser¬ 
vations at lower luminosities or with small exposures we grouped 
closely spaced observations with approximately the same averaged 
source intensity into three groups with the averaged luminosities of 
L 37 3.2, 6.2 and 12.7. The list of analysed observations with cor¬ 

responding averaged luminosities in the 3-100 keV energy range is 
presented in Table [T] These luminosity values may be regarded as 
close to the bolometric ones under the assumption that the bulk of 
the energy is released in the X-rays. 


HEXTE instruments onboard RXTE we used standard programs 
of the ftools/lheasoft v. 6.14 package and methods described in 
the RXTE cookboolQ. Depending on the spectrum quality photons 
were folded into 16, 8 or 4 phase bins, using individually deter¬ 
mined (per observation) and barycentrically corrected pulse peri¬ 
ods. 

The obtained spectra were approximated by a powerlaw model 
with an exponential cutoff (cutoffpl in the xspfc package) 

A£-^exp(-£/£cu,), (D 


where T is the photon index, £cut is the cutoff energy, A is a nor¬ 
malization. The spectrum was modified by one or several cyclolron 
absorption lines in the form 


exp 


-Tiv(£/£v)V^ \ 

(E-EMf + a-l)' 


( 2 ) 


To extract pulse phase-resolved spectra from the PCA and 


^ http://heasarc.nasa.gov/docs/xte/recipes/cook_book. html 
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Figure 2. Dependence of spectral parameters of V 0332+53 on the pulse phase for luminosities of L 37 = 36.9 (obs. 1, left) and L 37 = 12.7 (obs. 17, right). 
Corresponding pulse profiles are shown by black (upper panels) and grey (other panels) lines. Uncertainties correspond to the Icr confidence level. Two periods 
are shown for claiity. 


where En, <Tn and Tjv are the energy, width and depth of the cy¬ 
clotron line (N = 0) and harmonics (N > 0) , respectively (cy- 
CLABS in the xspec package; iMihara et alJI 199011 . The same model 
was used earlier for the analysis of pulse-phase averaged spectra 
and reasons for the choice of t he continuum and line models were 
discussed bv lTsygankov et aPlTOOdflOlOh . Note, that inclusion of 
the photoelectric absorption in the model does not affect the results 
of fitting in the 3-100 keV energy range and thus it was not taken 
into account. 

The first harmonic of the cyclotron absorption line was not 
confidently observed and its parameters were not well constrained 
in all data sets. Therefore, we examined possible effects of the in¬ 


clusion of higher harmonics in the model in such cases on the pa¬ 
rameters of the cyclotron line and of the continuum. Similarly t o 
the analysis of the pulse-averaged spectra dTsygankov et alJboObtl . 
the energy of the first harmonic Ei was assumed either to be a free 
parameter or fixed equal to the doubled energy of the fundamental 
Eq. We have found that the parameters of the continuum within the 
uncertainties do not depend on whether the higher harmonics are 
included in the model. The energy and width of the cyclotron line 
were found to slightly vary with the systematic uncertainty of about 
0.1 and 0.2 keV, respectively. In the following analysis these val¬ 
ues were taken into account when calculating uncertainties of these 
parameters. Finally, to take into account the iron fluorescent line at 
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~ 6.4 keV a Gaussian line at this energy with a fixed width 0.1 keV 
was included in the model as well. 

The total number of spectra which were analysed in the pa¬ 
per exceeds two hundreds and, therefore, it is not possible to fully 
discuss the results of fitting of the individual spectra here. We thus 
present in Fig. [T] several examples of pulse-phase resolved spectra 
for two luminosity levels to give some impression on the quality 
of the data and typical variation of spectral parameters with pulse 
phase. The best-fitting parameters are summarised in Table|2] It can 
be seen from Fig. [T] that the model adeq uately describes observa¬ 
tion s in a wide energy ran ge 3-100 keV. IPottschmidt et al.l l l2005ll 
and iNakaiima et aP J201(]|) have discussed a complex structure of 
the V 03324-53 spectrum near the cyclotron line. This is also vis¬ 
ible in the residuals of our spectra. These are likely related to the 
oversimplified description of the CRSF (i.e. with a Gaussian or a 
Lorentzian absorption profiles), whereas in reality it has probably 
a m ore complex shape (as discussed e.g. bv ISchdnherr et al]|2007l 
and iMukheriee & Bhattacharvall2012h . 


3 RESULTS 

The analysis had been carried out uniformly for all observations 
and a strong variation of both the continuum and the cyclotron line 
parameters with the pulse phase had been found at all luminosities. 
An example of such variability for two luminosities (L 37 = 36.9 
and L 37 = 12.7) is presented in Fig. together with the source 
pulse profile in the 3-20 keV energy band. To align the phases of 
individual observations we use the pulse profiles in 3-20 keV en¬ 
ergy range. Due to significant variation of the spin period and the 
pulse profile shape with luminosity, finding a robust timing solu¬ 
tion coherently defining pu lse phase at all luminosities is problem¬ 
atic dRaichur & Paulil2010h . Therefore, we roughly aligned pulse 
profiles obtained for individual observations to have the zero phase 
coincide with the main minimum of the pulse, and further on only 
discuss variations of the spectral parameters over the pulse profile. 
This approach does not require precise phase connection over the 
whole outburst. 

As can be seen from Fig. [2l the spectrum of the source is 
strongly variable with pulse phase. Moreover, the pulse-phase de¬ 
pendence of most parameters appears to be different at different 
luminosities (see detailed descriptions in following sections). We 
note that our spectral model has many parameters, some of which 
can be correlated (e.g., photon index versus cutoff energy, cut¬ 
off energy versus cyclotron li ne energy, see ICobum et alj 1 20021 : 
iLutovinov & Tsvgankovll201ll for results and discussion). There¬ 
fore, a change in one parameter might in principle be compensated 
by a change in another parameter with no change in the overall 
spectral shape. To demonstrate that this is not generally the case, 
and that the spectrum does indeed significantly change its shape, 
we present (as an example) in Fig.[3two spectra obtained at differ¬ 
ent pulse phases as well as their ratio. The variability both in the 
continuum and near the cyclotron line energy is apparent. Below, 
we discuss in details these variations. 

3.1 Continuum 

Currently there are no commonly accepted theoretical models 
to describe the spectra of accreting pulsars in details provided 
by observations, so empirical models are typically used instead. 
As it was noted above, the broadband spectra of V 0332-1-53 
(both the pulse-averaged and the phase-resolved) were described 



5 10 50 

Energy, keV 


Figure 3. (a) Energy spectra of V 03324-53 obtained for the 4th (i.e. phase 
0.1875-0.25; red points) and 7th bins (i.e. phase 0.375-0.4375; blue crosses) 
of the pulse profile (see Fig. [5) in the high luminosity state (L 37 = 36.9, obs. 
1, Table 1), and (b) their ratio. The solid lines represent best-fitting models 
for both spectra. 

with an absorbed power law with exponential cutoff at high 
energies. The study of the pulse-averaged spectra showed that 
the spectrum becomes softer and the cutoff energ y increases 
when luminosity drops jLutovinov & Tsvgankovl201lh . Moreover, 
the cutoff and cyclotron line energies were found to correlate 
jUntovinov & TsygankovluOl ih . although the physical reasons be¬ 
sides this correlation are not fully clear. In this paper, we use the 
same model as it is the simplest model which seems to capture all 
relevant features of the spectrum. 

3.2 Variations of the cyclotron fundamental line parameters 
with phase and luminosity 

As is evident from Fig.|2] the cyclotron line parameters signif¬ 
icantly vary with the pulse phase and the shape of these varia¬ 
tions (below we refer to them as ‘cyclotron energy/width/depth 
profiles’ for simplicity) depends on the source luminosity. The ob¬ 
served changes of these profiles should directly reflect changes 
in the structure of emission regions and their beam patterns with 
luminosity. To trace this variability we constructed the cyclotron 
energy/width/depth profiles for all observations from Table[T]and 
compared them with the corresponding pulse profiles. In Fig.|4l the 
pulse-phase dependencies of the cyclotron line parameters are pre¬ 
sented for eight levels of the source luminosity L37=36.9 (a), 32.8 
(b), 20.2 (c), 16.3 (d), 12.7 (e), 9.7 (f), 6.2 (g), and 1.7 (h). 

At high luminosities L 37 =^32-37 the cyclotron energy profile 
has a double-peaked shape with the maxima roughly coinciding 
with the maxima of the pulse profile. With decreasing of the lumi¬ 
nosity to L 37 =: 30 the profile becomes single peaked. Such a shape 
remains then more or less stable down to luminosity of L 37 =: 14. 
At even lower luminosities, the maximum of the cyclotron energy 
profile moves to the inter-peak minimum of the pulse profile and 
can be clearly seen for luminosities below L 37 10 . 

The evolution of the cyclotron depth profile with the lumi¬ 
nosity differs from the described above behaviour of the cyclotron 
energy profile. In particular, practically for all luminosities its max¬ 
imum roughly corresponds to the inter-peak minimum of the pulse 
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Figure 4. Variations of the cyclotron line energy (top left), width (top right) and depth (bottom) with the pulse phase for dilferent luminositie s: L 37 = 36.9 (a), 
32.8 (b), 20.2 (c), 16.3 (d), 12.7 (e), 9.7 (f), 6.2 (g), 1.7 (h) (black histograms). The relative columns heights according the reflection model jPoutanen et alJ 
l2Q13h are: hlR=0.91 (a), 0.84 (b), 0.57 (c), 0.48 (d), 0.4 (e), 0.32 (f), 0.22 (g), and 0.07 (h). Grey lines represent the corresponding normalized pulse profiles 
(right axis). The line width was fixed at 6 keV for the obseiwation with the lowest luminosity L 37 = 1.7 (obs. 23). 
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Figure 5. Dependence of the first harmonic energy (left) and the ratio of the line energies E\ /Eq (right) on the pulse phase for different luminosities: L 37 =36.9 
(a), 18.6 (b), 15.3 (c), 12.7 (d), 9.7 (e), 6.2(f) (black histograms). Grey lines represent corresponding pulse profiles. 


profile, and only at the very low luminosities covered by observa¬ 
tions (L 37 = 1.7) the maximum of the cyclotron depth profile shifts 
to one of pulse profile peaks. 

The cyclotron line width tends to anti-correlate with the flux 
at high luminosities (Fig.lUi, b, d) and to correlate with it at low 
luminosities (Fig.|4^, f, g). Such a behaviour is indeed expected in 
the reflection model. Generally speaking, although the width de¬ 
pends on several factors, it is expected to be broader in the case 
when we see the illuminated neutron star surface from aside and 
narrower when the illuminated surface is seen from the poles (see 
an extended discussion in Section 4.2). 


3.3 Variations of the energy of the first harmonic 

iTsygankov et al.l jlOOfih showed that the centroid energy Ei of the 
first harmonic roughly follows the luminosity dependence of the 
cyclotron line £o- Flowever, the ratio of the energies of the first har¬ 
monic to the fundamental one does change slightly with luminos- 
ity, from about ~ 2.2 at highest luminosities to ~ 1.9 at the lowest 
dTsygankov et al.l2006l : lNakaiima et al.l2010l) . Taking into account 
that during the pulse different parts of the emission regions are ob¬ 
served at different viewing angles and that the cyclotron line and 
the first harmonic could potentially be formed in slightly ditferent 
regions, it would be interesting and important to trace the behaviour 
of the first harmonic energy Ei and of the ratio £, /£o both with the 
pulse phase and luminosity. 

Note that the parameters of the first harmonic are reasonably 
constrained only in a subset of all obser vations even for pulse - 
averaged spectra (see detailed discussion in lTsygankov et al.l2006t) . 
The situation is even worse for the phase resolved analysis due 
to lower number of photons. Nevertheless, we found a number of 


observations, covering a luminosity range from L 37 =: 5 to =: 37, 
where a combination of the source intensity and exposure allowed 
us to determine and constrain the parameters of the first harmonic. 
Profiles of its centroid energy and corresponding ratio EJEq as a 
function of the pulse phase are presented in Fig.O 

It is difficult to make a detailed comparison of the Ei be¬ 
haviour with Eq, however, it can be seen that their profiles are not 
coincident (see Figs.|4]and|5]l. This is probably due to the fact that 
different parts of emission regions are responsible for the forma¬ 
tion of the cyclotron line and first harmonic of the CRSF. The ratio 
EijEo is also strongly variable (up to ~ 10 %) during the pulse. 
On the average, it demonstrates also a decrease with the luminos¬ 
ity, similar that found for p ulse-average spectra dTsygankov et al.l 
l2006l:lNakaiima et al.ll201ol) . 


4 DISCUSSION 


We now will try to describe qualitatively the observed pulse-phase 
changes of the cyclotro n line parameters usin g the recently pro¬ 
posed reflection model dPoutanen et al.l 20131). Acc ording to the 
widely accepted paradigm ( Basko & Sunyaevll 19761) . bright X-ray 
pulsars have relatively tall (comparable with the neutron star ra¬ 
dius) accretion columns at both magnetic poles. The columns are 
supported by the radiative pressure of the X-rays trapped within the 
optically thick accretion flow and slowly diffusing through the sides 
of the columns, so the plasma is most etfectively slowed down in 
the core of the column, where as the outer shell falls almost freely 
dUyubarskii & Svunvaev|[l983) . As a result, the emerging emission 
is expected to be strongly beamed towards the neu tron star surface 
due to the scattering on fast free-fal ling electrons dKaminker et al.l 
ll976l : lLvubarskii & Svunvaevlll988h . The formation of a cyclotron 
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a b 

Figure 6. Sketch of the angular distribution of different components for 
high (a) and low (b) luminosities. Red diagram corresponds to the direct 
component, blue one - to the reflected one. Only one accretion column is 
drawn for clarity. 


feature in the directly observed column spectrum is rather prob¬ 
lematic due to a strong gradient of the magnetic field along the 
column and scat t ering on free-falling electrons (see discussion in 
lNishimu^l2008l . l2014h . On the other hand, such a feature can be 
formed in the spectrum reflected from a n eutron star atmospher e 
where the magnetic field change is smaller l lPoutanen et al.ll201^ . 
The size of the spot on the neutron star surface illuminated by the 
accretion column depends on the column height, which, in turn, de¬ 
pends on the X-ray pulsar luminosity. The fraction of the column 
emission which is intercepted by the neutron star remains, however, 
large at all luminosities, so one does expect that a large fraction of 
the observed emission is reflected off the surface of the neutron star. 
It is not unlikely, therefore, that the observed cyclotron features also 
form in the atmosphere of the neutron star. This scenario allows to 
explain many important observational features of the CRSFs that 
otherwise are difficult to explain. For instance, the luminosity of 
V 0332-b53 in the considered outburst (and, therefore, the height 
of the accretion columns) changes by factor of twenty, whereas 
the line energy varies by just ~ 30 per cent, which is hard to ex¬ 
plain if the line is formed within the column. On the other hand, 
the dipole magnetic field over the neutron star surface decreases 
only by factor of two from the magnetic pole to the equator, so 
smaller variations of the CRS F energy are easier to accommodate. 
In fact, IPoutanen et al.l (12013h showed that the expected variation 
of the CRSF energy in this case is compatible with observations. 
Moreover, the reflection model was also capable of explaining the 
observed correlation of the line energy with luminosity. Below we 
will discuss qualitatively how the observed variation of CRSF pa¬ 
rameters with pulse phase and luminosity can be qualitatively ex¬ 
plained within this model. 


4.1 Phase dependence of the cyclotron line energy as 
function of luminosity 

The quantitative comparison of results of the phase-resolved spec¬ 
troscopy with the predictions of any theoretical models, including 
the reflection model, will be possible only after a self-consistent 
model emerges that describes the accretion column structure and 
the radiation field around the pulsar including the reflection from 



Figure 7. Dependence of the phase shift between maxima of the pulse pro¬ 
file and the cyclotron line energy on the X-ray pulsar luminosity. The 
error bars correspond to phase bins. 

the neutron star surface. Nevertheless some observational features 
described above could be understood qualitatively. 

If reflection off the neutron star surface is indeed responsible 
for the formation of the CRSF, the observed line parameters such 
as its energy, width and depth will depend on the size of illumi¬ 
nated spot(s). This, in turn, depends on many factors including the 
accretion columns height, beaming of the emerging radiation, the 
neutron star compactness, etc. Moreover, the observable part of the 
illuminated hot spots also depends on the relative position of the 
observer, being a function of the pulse phase. Thus, the cyclotron 
line parameters should be variable with the pulse phase. 

From geometric considerations, the line energy is expected to 
be maximal when one of magnetic poles is directed to the observer, 
so that the polar regions with the strongest field are observed. At 
other angles large fraction of reflected radiation will be scattered 
further from the poles, where the field is lower, and a slight decrease 
in the observed CRSF energy is expected. 

It is interesting to note that the observed line energy maxi¬ 
mum tends to coincide with one of the total flux maxima at high 
luminosities (see Figs |4^-e) and it tends to be between two flux 
maxima at low luminosities (Figs|4f-h). A quantitative analysis this 
behaviour is not simple. Pulse profiles as well as the dependence 
of the line energy on the phase are very different and variable. 
This makes the standard correlation analysis problematic. Instead, 
a phase shift between the maxima of the considered curves can be 
used to quantify the potential correlation. To account for the appar¬ 
ent asymmetry of the profiles, we used the centre of the segment 
connecting rising and declining parts of the pulsed maximum at 
half of the maximum flux as a reference phase. Results of such an 
analysis are presented in Fig.|7]and they are in agreement with our 
expectations - the phase shift is a relatively small (< 0.1) at high 
luminosities and increases at low luminosities. 

In the framework of the reflection model, the cyclotron line 
only contributes to the reflected spectrum, whereas total flux at 
a given pulse phase includes both reflected and directly observed 
emission. Therefore, the observed change in the pulse phase be¬ 
haviour of the cyclotron line energy can be explained if the ratio 
between these two components changes with luminosity. The ratio 
of the reflected emission to that observed directly from the column 
depends on total fraction of intercepted emission and on the beam¬ 
ing characteristics of the two components. The reflected compo¬ 
nent is expected to have a broad pencil-like beam pattern directed 


























Pulse phase spectroscopy of V 0332+53 9 


along the magnetic dipole axis. The beam pattern of the emission 
observed directly is more complicated and depends on luminosity. 

At high luminosities, both columns are expected to contribute 
to the observed flux at all pulse phases with the combined beam pat¬ 
tern similar to that of the reflected component (see Fig.[^). There¬ 
fore, the pulse profile is expected to have two maxima correspond¬ 
ing to the two illuminated spots around the two magnetic poles. In 
such a situation, the maximum of the cyclotron line energy profile 
will coincide with one of the pulse profile maxima, which in turn 
corresponds to the pole with the smallest inclination angle to the 
line of sight. 

At lower luminosities, the height of the column decreases, and 
the illuminated spot around the second pole becomes less visible 
or even fully obscured by the neutron star. As a result, a fan-like 
beam pattern from the two columns is expected in this case to have 
two maxima when both columns are observed sideways (see Fig.[^. 
The reflected component retains the pencil-like beaming and, there¬ 
fore, is expected to be in anti-phase with the direct emission. Even 
though the total fraction of the reflected emission is expected to in¬ 
crease at lower luminosities, only contribution from one pole can 
be observed at a time. On the other hand, both columns are ex¬ 
pected to be visible when the contribution of the directly observed 
component is maximal, so the overall maximum of the pulse profile 
shall be defined by the component observed directly. Still, the cy¬ 
clotron line is only present in the reflected component, which is in 
anti-phase with the directly observed emission. Therefore, the max¬ 
imum of the cyclotron line energy profile is expected to move away 
from the pulse-profile maximum at lower luminosities in agreement 
with observations. 

We can find a relative column heigh t in the framework of 
the reflection model fPoutanen et alJl^lSh . We used equation (2) 
in that paper and substituted model parameters which give the 
best fit to the observ ed dependence ii(,yc versus L (see Fig. 5 in 
IPoutanen et al.ll 201 ^ . We find that the relative column height hJR 
decreases from 0.89 for the brightest state to 0.07 for the dimmest 
one (the relative column heights for all the luminosities are given 
in the caption to Fig.|4ll. In framework of this model, the first de¬ 
scribed scenario with the pencil-like reflection component domina¬ 
tion is realised during the high-luminosity states and relatively tall 
columns with hIR > 0.5. The second scenario, with the direct fan¬ 
like component, is likely realised for short columns with hIR < 0 . 2 . 

Finally, we note that the light bending is an im portant factor 
influencing the visibility of illum inated spots (see e.g. lseloborodovl 
l2002t lAnnala & Poutan^ll999h which we did not account in out 
qualitative discussio n. However, estimation s of h/R are done on the 
base of the model bv IPoutanen et alj(l2013h who accounted for light 
bending. 


4.2 Correlations between the parameters of the cyclotron line 

Based on t he _ ^ectral analysis of several X-ray pulsars, 

ICobum et al.l (l 2002 h found that there are two types of correla¬ 
tions between cyclotron line parameters: between its width cto and 
energy and between the relative cyclotron line width (TqIEo 
and the line depth Tq. Late r, similar correlations were reported by 
iKrevkenbohm et ^ ( |2004) for pulse phase-resolved spectra of the 
X-ray pulsar GX 301-2. At the same t ime no obvious correlations 
were found bv lTsygankov et afl ( 1201 Oh from the analysis of pulse- 
average spectra of V 0332-1-53 in a wide energy band. 

To obtain the relations ctq - Eq and (ctoIEo) — To we have 
reanalysed uniformly all data sets from Table [T] dividing the cor¬ 
responding pulse profiles into eight phase bins and approximating 


the obtained spectra with the model described above (cutoffpl + 
CYCLABs). The observations where cyclotron line parameters were 
not reasonably constrained due to insufficient statistics (effectively 
limiting the number of phase bins to four or less) were excluded 
from further analysis. Obtained results are presented in Fig.[^ In 
general, its left panel is similar to the figure 5 from the paper of 
iTsygankov et al.l l l2006h . where the dependence of the cyclotron 
line width on its energy was shown for pulse phase averaged spec¬ 
tra. Note a fairly strong scatter of the data points for large cyclotron 
energies, which is probably caused by lower statistics at low lumi¬ 
nosities. 

Note that we used the cyclabs model to describe the cy- 
clotron absorption li ne as opposed to m a ny other authors (e.g., 
Coburn et~^ l2002l : IKrevkenbohm et af] 1 20041 : iKlochkov et al.l 
20081 ) . who used the gaussian absorption line profile (gabs in the 
xsPEC package): 



(3) 


where Ei^, ctn, and Tn are the energy, width and optical depth of 
the cyclotron line (N = 0) and harmonics (N > 0), respectively. 
Both models adequately approximate cyclotron absorption line and 
its harmonics (see Table[3] where best-fitting parameters for both 
CYCLABS and gabs models are presented). Note, that the cyclotron 
line energy derived from the cyclabs model is systematically lower 
(by ~l-3 keV) than the energy derived from the gabs model ( see 
e.g. lMih^ll995l : [Nakaiima et alJl201^ : lTsvgankov et alj2012l) . 

In this paper, we attempt not only to find spectral parameters 
of V 0332-1-53, but also to search for their dependence on the pulse 
phase and luminosity and possible physical correlations between 
them. Therefore, a correlation analysis was carried out to investi¬ 
gate potential intrinsic correlations of CRSF parameters for cyclabs 
and GABS models. To illustrate obtained results, we plot in Fig.j^the 
confidence contours of the cyclotron line width <To versus cyclotron 
line energy Eq obtained with cyclabs and gabs models for luminos¬ 
ity L 37 =; 15.3 (obs. 15) (results for other luminosities are similar). 

Several conclusions can be made from this analysis: the line 
energy depends on the CRSF models (see above); an accuracy of 
the line energy measurements is better for the cyclabs model, while 
the accuracy of the line width measurements is better for the gabs 
model; there is an intrinsic correlation between ctq and Eq for the 
GABS model, while no correlation and fairly significant scatter of the 
data points or even anti-correlation are seen for the cyclabs model. 

A linear regression fit to the data shown in both panels of 
Fig.[8]gives a hint that there is indeed some weak anti-correlation: 
0-0 = (12.7 ± 1.6) - (0.24 ± 0.06)£o and cto/Eo = (0.30 ± 0.03) - 
(0.030 ± 0.015)t o . Thu s, we confirm here the previous result by 
ITsygankov et alj ( l2010h that there is no positive correlation be¬ 
tween the cyclotron line parameters for V 03324-53 if the cyclabs 
model is used for the cyclotron line. Finally, we note that if the 
faintest observation is excluded from Fig.[^ (black squares), the 
significance of the anti-correlation is increased and in this case 
0-0 = (17.0 ± 1.2) - (0.41 ± 0.05)£o and cto/Eq = (0.42 ± 0.03) - 
(0.095 ±0.016)to. 

Formation of the CRSF in the framework of the reflection 
model is still pretty much work in progress, so here we restrict 
ourselves to qualitative discussion of the results. The cyclotron line 
width ctq is mainly determined by two factors. First of all, the line is 
thermally broadened by a factor which depends on the plasma tem¬ 
perature in the line forming region and the angle between photon 
momentum direction and the magnetic field. The maximum ther¬ 
mal line width is realised for the radiation propagating along the 
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Figure 8. (a) Cyclotron line width (Tq versus centroid energy Eq,. (b) Relative cyclotron line width <to/-E'o versus its depth tq. Different colours and symbols 
correspond to observations with different luminosities (see Table^: obs. 1 — red squares, 4 - black triangles, 5 - green circles, 8 - blue squares, 13 - magenta 
circles, 15 - green triangles, 17 - cyan circles, 18 - black circles, 19 - cyan triangles, 20 - red circles, 21 - black squares. 




Figure 9. Cyclotron line width (Tq vs. cyclotron line energy E(^, obtained using spectral fits with cyclabs (left panel) and gabs (right panel) models for 
phase-resolved spectra of V 0332+53 in observation 15 (Table^. Contours (solid, dashed, dotted lines) correspond to the 1,2,3 cr levels, respectively. 


Table 3. Pai'ameters of the best fits with cyclabs and gabs models to the same pulse phase-resolved spectra of V 0332+53 (Obs. 15, L 37 = 15.3). 


Model 

r 

Fcut 

Eq 

O-Q 

tq 

El 

o-\ 

r\ 

(d.o.f.) 



keV 

keV 

keV 


keV 

keV 







1st bin (phase 0.0-0.125) 





GABS 

0.38 ± 0.02 

7 4Q+0.20 

z,:?.uz,_o 14 

c o 1 +0.18 
-0.08 

1 67+0.07 
-0.03 

52 2+®-’^ 
- 0.2 

- 0.2 

2.45 ± 0.11 

105.9(107) 

CYCLABS 

0.32 ± 0.03 

8.36!“3 

26.65 ± 0.10 

o.ju_o 12 

1.73 ±0.04 

51.6 ± 0.6 

^■^- 0.7 

^■^-0.7 

109.2(107) 





7th bin (phase 0.75-0.875) 





GABS 

0.36 ± 0.04 

8.56 ±0.15 

27.90 ±0.15 

5.04 ±0.09 

1.93 ±0.07 

54.8 ± 1.0 

10.5 ± 0.2 

2.8 ± 0.3 

108.3(107) 

CYCLABS 

0.30 ± 0.04 

10.14 ±0.25 

25.85 ±0.08 

^“^^-0.06 

2.09 ± 0.04 

52.2 ±0.7 

10.5 ± 1.1 

3.2 ± 0.3 

107.2(107) 
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Bolomelric luminosity, 10®'’’ erg s"‘ 


Figure 10. (a) Maximal difference between cyclotron line energies over the 
pulse versus source luminosity, (b) Pulse fraction in the energy range 25- 
45 keV versus bolometric luminosity for V 0332+53 (from Tsygankov et 
al., 2010). 


magnetic field A£d ~ E (2kT or about 4 keV for the X- 
ray pulsar V 0332+53 (E « 25 keV, kT « - 5-10 keV). It 

is also well known and might be important for the future quanti¬ 
tative analysis, that the cyclotron line op acity strongly depends o n 
the direction as well (see e.g. Fig. 1 in ISuleimanov et alj|2012h . 
The second important factor is a dispersion of the magnetic field 
strength over the line forming region. We focusing helow on this 
factor. According to the reflection model, the CRSF originates in 
radiation reflected by the NS atmosphere at different latitudes. The 
magnetic field decreases from poles to equator by factor of two, so 
the dispersion of the field within the line forming region depends 
on the size of the illuminated hot spot around the magnetic pole. 
If we assume that the thermal broadening does not depend on the 
accretion column height, the line width will increases with the pul¬ 
sar lumi nosity, while the cyclo tron line energy will decrease (see 
Fig. 4 in IPoutanen et al.ll2013h . Indeed, this agrees with observa¬ 
tion (see Figs|4]and[8ji). Strong deviations from these dependences 
take place only at low luminosities when the relative column height 
is quite small hjR < 0.2 (the black squares in Fig.[^. In this case 
the situation can be more complicated because the cyclotron line 
can be formed not only on the stellar surface but also in the small 
accretion column and/or (prohahly hotter) plasma deceleration re¬ 
gion. 


4.3 Variations of the cyclotron line energy dnring the pulse 
as fnnction of luminosity 

To character quantitatively the variability of the cyclotron line en¬ 
ergy with the pulse phase we introduce a parameter AFo ^ax, which 
is the difference between maximum and minimum values of the cy¬ 
clotron line energy Eq over the pulse. If we expect that during the 
pulse different parts of the emission regions are observed, this pa¬ 
rameter should be related to the geometrical size of these regions 
and should reflect their changes with variations of the source lu¬ 


minosity. In the case of super-critical sources (with the accretion 
column above the surface), the parameter AFo.max is related to the 
column height and size o f the illuminated part of the neutron star 
surface, correspondingly jPoutanen etalj2013h . 

The dependence of AFo.max on the luminosity is shown 
in Fie.llOh. Despite of some scattering, we see that AFo^max 
decreases with increase of the luminosity up to L37 « 30. At higher 
luminosities Aito.max increases quickly practically up to values, 
typical for low luminosities. To quantify this, we approximated the 
A£^o,raax “ E^-j dependence with the linear fit in different luminosity 
ranges. We found that AFo.max = (1.65 ± 0.20) - (0.03 ± 0.01)L37 
for L 37 < 30 and AFo.max = (-2.4 ± 1.8) + (0.10 ± 0.05)L37 for 
^37 > 30. This gives > 99 per cent significance that the difference 
in slopes is real. Finally, we note that although it is hard to interpret 
such behaviour of AiTo.max, it is similar to the behaviour of the 
pulse fraction with the source luminosity (Fig. QSt). 


5 SUMMARY 

The pulse phase-resolved spectroscopy of X-ray pulsars provides a 
lot of additional information in a comparison with the analysis of 
the pulse-averaged spectra, therefore it can be used as a powerful 
instrument for the discrimination between different models for the 
X-ray pulsars emission mechanisms. Up to date the pulse phase- 
resolved spectroscopy has been done for several sources. However, 
until now the limitations of observational data did not allow to in¬ 
vestigate the influence of the source luminosity which determines 
geometrical structure of the emitting region. 

In this paper we presented results of the pulse phase-resolved 
spectroscopy of the transient X-ray pulsar V 0332+53 based on the 
RXTE observations carried out during the powerful outburst in Dec 
2004 - Feb 2005. For the first time such a study was performed for 
the same source in a wide luminosity range L 37 (1 - 40). 

It was shown that both continuum and cyclotron line param¬ 
eters are strongly variable with the pulse phase at all luminosities. 
The most important conclusion is that all spectral parameters are 
variable not only with the pulse phase, but the pattern of this vari¬ 
ability with respect to the pulse profile is different for different lu¬ 
minosities. Such a behaviour reflects directly changes in the struc¬ 
ture of emission regions and their beam pattern with the luminosity 
changes. 

For instance, at very high luminosities L 37 (32-37), the cy¬ 

clotron energy profile has double-peaked shape almost coinciding 
with the pulse profile. With the decrease of the source intensity to 
L 37 30 the cyclotron energy profile is transformed to the single 

broad peak with the maximum at the same phases where the pulse 
profile peaks. Such a shape remains more or less stable till the lumi¬ 
nosity of L37 =: 14. At even lower luminosities below L37 10 the 

maximum of the cyclotron energy profile moves to the inter-peak 
minimum of the pulse profile. 

We interpreted these observational results qualitatively in the 
framework of the recen tly proposed reflectio n model for the cy¬ 
clotron line formation l lPoutanen et alj|2013h . Nevertheless, it is 
necessary to say that deep observations with more sensitive instru¬ 
ments in hard X-rays (e.g., NuSTAR) are needed to make quantita¬ 
tive concl usions and verify an applicability of competitive models 
(see, e.g.. lNishimurall20141 . 
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